Single Stranded DNA Translocation Through A Nanopore: 
A Master Equation Approach 
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We study voltage driven translocation of a single stranded (ss) DNA through a membrane channel. 
Our model, based on a master equation (ME) approach, investigates the probability density function 
[pdf] of the translocation times, and shows that it can be either double or mono-peaked, depending 
on the system parameters. We show that the most probable translocation time is proportional to 
the polymer length, and inversely proportional to the first or second power of the voltage, depending 
on the initial conditions. The model recovers experimental observations on hetro-polymers when 
using their properties inside the pore, such as stiffness and polymer-pore interaction. 
PACS number: 87.14.Gg, 87.15.Aa,87.15.He 



1. INTRODUCTION 

Translocation of biopolymers through a pore embed- 
ded in a membrane is a fundamental step in a variety 
of biological processes. Among the examples are the 
translocation of some species of m-RNA through the nu- 
cleus membrane which is the first stage of gene expres- 
sion in eucaryotic cells [1], and the attack of cells by 
viruses that occurs by injecting the genetic information 
through a hole in the cell membrane [2]. The translo- 
cation importance in bio-systems, and the possibility for 
developing fast sequencing methods, have been the moti- 
vation for recent experiments, in which a voltage-driven 
ssDNA is translocated through a-hemolysin channel of 
a known structure [3]. In these experiments one mea- 
sures the time it takes a single ssDNA molecule to pass 
through a membrane channel [4-7]. Since ssDNA is neg- 
atively charged (each monomer has an effective charge of 
zq, where 2; (0 < z < 1 ) is controlled by the solution 
pH and strength), when applying a voltage the polymer 
is subject to a driving force while passing through the 
membrane from the negative (cis) side to the positive 
(trans) side. Because the presence of the ssDNA in the 
transmembrane pore part {TPP) blocks the cross- TPP 
current, one can deduce the translocation times prf/from 
the current blockade duration times [4-7]. It has been 
found that the shape of the translocation times pdf is 
controlled not only by the voltage applied to the system, 
the temperature and the polymer length but also by the 
nature of the homo-polymer used: poly-dA (A- adenine); 
poly-dC (C-cytosine); poly-T„u (r„„-thymine), and the 
sequence of hetro-polymers [4-7]. 

The translocation process can be roughly divided into 
two stages. The first stage is the arrival of the polymer 
in the vicinity of the pore and the second stage is the 
translocation itself. Several models have been suggested 
for describing the translocation stage. In [8] an equation 
for the free energy of the translocation, obtained from 
the partition functions of the polymer parts outside the 



TPP, was derived and used to calculate the mean first 
passage time (MFPT). Other investigators used similar 
ideas with improved free energy terms by taking into con- 
sideration effects such as the membrane width [9-10], or 
assumed that only the part of the ssDNA inside the TPP 
affects the dynamics of the translocation rather than the 
polymer parts outside the TPP [11-12]. 

In this work we present a new theoretical approach 
that allows to consider both the polymer parts outside 
the TPP and within the TPP. Using the ME we are able 
to map the three dimensional translocation onto a dis- 
crete space one-dimensional process. Based on the ME 
we compute the pdf of the first passage times (FPT) of the 
translocation, F{t), and the MFPT, as a function of the 
system parameters. We relate our theoretical results to 
recent experimental observations and by analyzing them 
using our model, we come up with physical understand- 
ing of these observations. 

2. THE MODEL 

An n(= N + d — l)-state ME is used to describe the 
translocation of an iV-monomer long ssDNA subject to 
an external voltage V, and temperature T, through a 
TPP of a length that corresponds to d{— 12) monomers. 
A state is defined by the number of monomers on each 
side of the membrane and within the TPP. A change in 
the state of the system (only nearest states transitions are 
allowed) is assumed to be controlled mainly by the be- 
havior of the polymer within the TPP in the presence of 
the applied voltage. In addition, it is assumed to be influ- 
enced by entropic and interaction factors of the polymer 
outside and within the TPP. Absorbing ends are used 
as boundary conditions, which are the natural choice for 
this system because the polymer can exit the TPP on 
both sides. The state j — n represents the arrival of 
the first monomer into the TPP from the cis side, and 
the state j = represents the departure of the far end 
monomer from the trans side of the TPP. Let Pj (t) be 
the pdf oi state j that contains rrij nucleotides occupying 
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the TPP according to: 

r J, j<d (1.1) 

mj = I d, N > j > d given N > d, (1.2) 
[n + l-j, j>N (1.3) 

and a similar set of equations for a short polymer, N < d, 
which is obtained by exchanging N and d in Eqs.(l.l) 
through (1.3). The governing equations of motion are: 

dPj{t)/dt = aj+ijPj+i{t) + a,_i,,P,_i(t) 

-{ajj+i + aj.j-i)Pj{t); .7 2, • • •, n - 1, 

dPy{t)/dt = ay+l^yPy+l{t)Sy 

V,n 

-{ay^y-l+ay^y + l)Py{t); ^=1,72. (2) 

Equations (2) can be written in matrix representation 
dP /dt = aP . The propagation matrix A, is a tridiag- 
onal matrix that contains information about the transi- 
tions between states in term of rate constants. We as- 
sume that the rate constants can be decoupled into two 
terms: 

^ ko{T)P3.j+/-iiV^T)- (3) 

The first term provides the rate to perform a step, while 
the second term gives the probability to jump from state 
J in a given direction, and obeys: Pj.j+i +Pj.j-i = 1. To 
obtain kj , we first consider the bulk relaxation time of a 
polymer [13] n oc P^bb'^N^ , where /^-^ = ksT, 6 is a 
monomer length, ^t, is the Stokes bulk friction constant 
per segment, = 677677 (where r] is the solvent viscosity), 
N is the number of monomers in the polymer and the 
dimensions dependant //, represents the effect of the mi- 
croscopic repulsion on the average polymer size. In three 
dimensions /U = 3, 9/5, 3/2 for rod-like, self-avoiding and 
Gaussian (Zimm model) chains, respectively. To com- 
pute the relaxation time inside the TPP, the confined 
volume of the TPP must be taken into account. For 
a rod-like polymer the restricted volume dictates a one- 
dimensional translocation, whereas for a flexible polymer 
the limitations are less restricted. We implement this by 
taking as a measure of the polymer stiffness inside the 
TPP that obeys: < /x < 1.5. The expression for the 
relaxation rate of state j is therefore, 

k, = l/{f3Cpb^m^) = R/m';. (4) 

From Eq.(4) it is clear that as /x becomes smaller the 
rate to preform a step becomes larger, namely kj for a 
rod-like polymer increases. This feature appears, at first 
sight, to be in contradiction to the relaxation time be- 
havior of a bulk polymer, where a rod-like polymer has 
a larger relaxation time than that of a flexible polymer. 
This contradiction is resolved by taking into account the 
diflferent dimensional demand for a rod-like polymer rel- 
ative to a flexible polymer inside the TPP. Because /i is 
a measure of the polymer stiffness inside the TPP, it is 



controlled by the interaction between the monomers oc- 
cupying the TPP, e.g base stacking and hydrogen bonds, 
and therefore is affected by the monomer type and the 
sequence of the ssDNA. 

The friction constant per segment inside the TPP, ^p, 
represents the interaction between the ssDNA and the 
TPP. The physical picture is that during translocation 
there are few or no water molecules between the polymer 
and the TPP. Consequently the water molecules inside 
the TPP can hardly be viewed as the conventional sol- 
vent and the Stokes friction constant is replaced by 
representing the ssDNA- TPP interaction. 

To calculate the second term on the right hand 

side of Eq.(3), we assume a quasi-equilibrium process and 
use the detailed balance condition for the ratio of the rate 
constants between neighboring states: ajj_i/aj_i,j = 
e-/3'^^j, where AEj = Ej_i - Ej. We then use the 
approximation ajj-i/aj-ij ~ pj — pjj-i), and 
deduce the jump probabilities by computing AEj. To 
compute Ej the contributions from three different sources 
are considered: electrostatic, Ej, entropic, Ej, and an 
averaged interaction energy between the ssDNA and the 
pore, 

For the calculations of the electrostatic energy differ- 
ence between adjacent states, ^Ej, we assume a linear 
drop of the voltage along the TPP and write for rrij 
monomers occupying the TPP penetrating from the cis 
side of the membrane for a length of ijb: 

ruj +ij —1 

EP^zq{V/d) ^ n = zq{V/2d)mj{mj +2ij -I). 

n=ij 

(5) 

The effective charge per monomer in the TPP is taken 
to be the same as of the bulk. For states that contain 
monomers at the trans side of the membrane, zqV should 
be added for Ai^J. This contribution represents the addi- 
tional effective charge that passed through the potential 
V. Accordingly, the expression for Ai^J is (see Appendix 
A): 

Ai;j = zqV{mj + aj)/d. (6) 

Here aj takes the values: aj = { — 1;0;1} for cases 
described by equations {(1-1); (1-2); (1-3)} respectively 
{aj = 1 and aj = —1 correspond to the entrance and 
exit states of the translocation, whereas aj = corre- 
sponds to the intermediate states of the translocation). 
For a short polymer, aj has the same values as for a large 
polymer. 

The contribution to I^Ej from A£^| is composed of 

two terms. One term is the entropic cost needed to store 
rrij monomers inside the TPP, and the second term orig- 
inates from the reduced number of conflgurations of a 
Gaussian polymer near an impermeable wall. Combin- 
ing these terms leads to (see Appendix B): 

Ai?; = 7,7/3, (7) 
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where 7^ = { — 1 + g^; 5^; 1 + g^}, for cases described by 
equations {(1-1); (1-2); (1-3)}, respectively, gj is given in 
Appendix B in terms of Nj^^is and Nj^trans^ which are 
the number of monomers on the cis and trans sides cor- 
respondingly. For a short polymer jj behaves similarly 
but for intermediate states gj=0. 

For computing Ai?* we focus on the average interac- 
tion between the ssDNA and pore (not only its trans- 
membrane part). Due to the asymmetry of the pore 
between the cis and the tran sides of the membrane 
[3], the ssDNA interacts with the pore on the cis side 
of the membrane and within the TPP but not on the 
trans side of the membrane. Assuming an attractive 
interaction, Ej = —ksT^N — Nj^trans) (more properly, 
Ej = —s{N — Nj^trans)^ and we set e = fcsT in the rele- 
vant temperature regime) we obtain: 



^E] = Cj/13, 



(8) 



where Q = {1; 1; 0}, for the cases described by equations 
{(1.1); (1.2); (1.3)}, respectively, and for a short polymer 
= {1;0;0}. 

Using the above relations, and defining Sj = -yj + Q we 
obtain 



= (1 + e' 



(9) 



For the system to be voltage driven P\AEj\ > Sj must be 
fulfilled, which translate into the condition: V/Vc > 1, 
where a characteristic voltage is introduced: Vq^ = 
(1 + l/d)l3z\q\. This inequality ensures that there is a 
bias towards the trans side of the membrane. Otherwise 
the polymer is more likely to exit from the same side it 
entered than to transverse the membrane. Under exper- 
imental conditions [6] Vc = 46rnF, when using z « 1/2. 

In Figure 1 we show the different contributions to AE'j: 
/3Ai;J (for Pz\q\V=l) and 5j. /SAE^ decreases for the 
entrance states of the translocation, increases at the exit 
states of the translocation, and is a negative constant 
for intermediate states. Clearly I3AE^- < reflects the 
field directionality. On the other hand, 6j opposes the 
translocation for the entrance and intermediate states. 
For the entrance states 5j > due to both entropic terms, 
but approaches zero (from below) for the exit states of the 
translocation, due to the cancellation of AEj against the 
entropic gain of storing less monomers within the TPP. 
At intermediate states 5j » 1, where its shape near the 
crossover between the different situations is controlled by 
9j- 

3. RESULTS AND DISCUSSION 

3.i. The FPT pdf 

In this subsection we compute the FPT pdf, F(t), and 
examine its behavior as a function of the system param- 
eters. F{t) is defined by: 
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FIG. 1: The components of the AEj, (3AE^^ {(3z\q\V = 1) 
and 5j , are shown for TV = 100. 



where the survival probability, namely, the probability to 
have at least one monomer in the TPP, S{t), is the sum 
over all states pdf 



S{t) = TfCe^^C-^Vo. 



(11) 



F{t) = d{i - s{t))/dt. 



(10) 



Here if is the summation row vector of n dimensions, 
is the initial condition column vector, (Pa)j = Sx.j, 
where x is the initial state, and the definite negative real 
part eigenvalue matrix, D, is obtained through the trans- 
formation: D = C~^AC. 

Substituting Eq.(ll) into Eq.(lO), F{t) is obtained. 
Figure 2 shows the double-peaked F(t) behavior as a 
function of V/Vc and N (inset), for starting at x = 
N + d/4:. The left peak represents the non- translocated 
events. Its amplitude decreases as V/Vc increases but 
remains unchanged (along with the position) with the in- 
crease in N, because only the 'head' of the polymer is in- 
volved in these events. As x decreases, namely when the 
translocation initial state shifts towards the trans side, 
the non-translocation peak and the translocation peak 
merge and for x = n/2, namely, for an initial condition 
for which Nj_cis=Nj,trans, F{t) has one peak (data not 
shown) independent of V/Vc- 

For a convenient comparison to experimental results, 
we calculate and examine the behavior of the most prob- 
able average translocation velocity, v = xb/tm, where 
is the time that maximizes the translocation peak. 

Figure 3 shows the v(N) behavior, for V/Vc = 
2.60. v{N) is a monotoniacUy decreasing function of N 
for N < d, but is independent of N for N > d. This 
agrees with the experimental results [6]. The origin for 
v{N) behavior, stems from Eqs.(l) and (4), according to 
which the minimal kj is a decreasing monotonic function 
of N for N < d, but is independent of N for N > d. 

The velocity, v{V), depends on x, v{V) = Vx{V), and 
changes from a linear to quadratic function of the voltage 
when changing x, as shown in the inset of figure 3, for 
N = 30. 
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FIG. 2: F{t), for several values oiV/Vc, with: N=30, x=N+ 
d/4, z « 1/2, T=2°C, Vc=46mV, //=1.14 and i?=10^Hz. 
The left peak represents the non-translocated events, whereas 
the right peak represents the translocation. Inset: F{t) as a 
function of A'', with V/Vc = 1.5, and the other parameters as 
above. The non-translocation peak remains unchanged when 
increasing the polymer length, while the translocation peak 
shifts to the right with N. 



Starting at x = A'' + 1, a linear scaling is obtained: 

vn+i{V) = bi{V — 62/^1)- The coefficient 62/'^! can be 

^ — 1 

identified as an effective characteristic voltage: Vc = 
(1 + l/d)Pz\q\ = h\/b2- From the last equality z can be 
extracted. 

Starting at x = N/2, i.e. when initial state is close 
to the exit states, a square dependence is obtained: 
vn/2{V) = ci{V - C2)2 + C3, with ci = o(10-5), C2 = 
40mV and C3 = o(10^^). These coefficients are similar 
to the measured values [6]. 

We note that both linear and square scaling behav- 
iors have been observed experimentally [4,6]. A possible 
explanation for the different functional behavior of v{V) 
might originate from different data analysis that can be 
interpreted as having a different initial condition. 

To get numerical values for and fx, we use 
the experimental data in refs. [6-7], and ob- 
tain: ^(C) = 1, fi{A) = 1.14, iiiTnu) = 1.28, and 
^p(A) « ^C) = Cp(T„„) = Cp(^)/3. 

Prom these values we find the limit in which the relax- 
ation time of the polymer parts outside the TPP can 
be neglected. We estimate the maximal bulk number 



N„ 



for which the bulk relaxation 



of monomers, 
time is much shorter (5%) than of the TPP relax- 
ation time. For a po\y-dA bulk Zimm chain we get 
Nmax ~ 271, by taking for the viscosity the value 
for water at 2°C, r; w 1.7 • lO-^Ns/rn^. Using this 
value, can be calculated from the Stokes relation 
to be: ~ 'i-0~^ meVs/nrn^, which is three order of 
magnitude smaller than ^p. 
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FIG. 3: The velocity as a function of the polymer length, for 
the same parameters as in figure 2, various initial conditions 
(linear-log scales) and V/Vc = 2.6. v{N) tends towards a 
length independent, but display a sharp decrease for short 
polymers. Note for the step feature of v{N) when N becomes 
larger than the TPP length, d. Inset: v{V) , for iV = 30, 
showing linear and quadratic scaling depending on the initial 
condition, x. 



3.ii. The MFPT 

Additional information about the translocation can be 
obtained by computing the MFPT, r. To compute an 

analytical expression for t, we consider a large polymer, 
N > d, and replace Pjj-i and kj by state independent 
terms: p+ = [1 + e(-^/^c+i)]-i, valid for x ^ N, and 
k = l/{/3£^pb'^d^). This leads to = p+k and a_ = 
(1— p_i_)fc, which defines a one dimensional state invariant 
random walk. The MFPT is obtained by inverting A 
[14]: r = J^tF{t)dt = -if A-'^Vq. The calculations 
of the elements of the state independent yields (see 
Appendix C): 



A(p")A(p"+i-")p^'" 
A(p)A(p"+i) k ' 



s <x, (12) 



where A(p"') = — and (— A"-^)^,^; for s > x is 
obtained when exchanging s with x and with p_ in 
Eq.(12). Summing the x column elements of (— A~^) we 
obtain r (see Appendix C): 

A(p"+i-^)p^x - A(p^)j3!::+i-^(n +l-x) 
^~ fcA(p)A(p"+i) ' 

which in the hmit of a weak bias, V/Vc ^ 1, can be 
rewritten as (see Appendix C) 



2x£,ph'^d'' 



z\q\{l + l/d)V -Vc 



(14) 



Although r and tm are different characteristics of F{t) 
and differ significantly when slow translocation events 
dominate, Eq.(14) captures the linear scaling with 
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TV and XjV . The quadratic scaling of tm with \jV 
is obtained when using Eq.(9) rather than its state 
invariant version, for starting at, or near, an initial state 
for which < 0. 

3.111. The sequence effect 

Under the assumptions that are presented below, we 
now construct and \x for every ssDNA sequence and 
thus examine the sequence effect on t„-i. For a given ss- 
DNA sequence occupying the TPP in the j state, we write 
an expression for the average friction of that state, ^pj, 
assuming additive contributions of the monomers inside 
the TPP: 



calculation of fxj following the definition: 



1 



(15) 



Here jiUs stands for the nucleotide ,s occupying the 
TPP. To construct a compatible state dependent stiff- 
ness parameter, /Zj, we first argue that only nearest 
monomers can interact inside; the TPP and thus con- 
tributes to the rigidity of the polymer, which in turn 
increases kj. We then examine the chemical structure 
of the nucleotides and look for 'hydrogen-like' bonds be- 
tween adjacent bases. The term 'hydrogen-like' bonds 
is used because the actual distance between the atoms 
that create the interaction may be larger than of a typ- 
ical hydrogen bond. The pairs AA and CC can interact 
but not the pairs TnuTnu and CA. For the pair CT„„ 
the interaction is orientation dependent; namely, for the 
I pairs sequence, poly-<i(CT„„);, the interaction is within 
each of the pairs but not between the pairs. Accordingly 
we have: pi{nu nu) = fJ,{nu), IJ,{CA) = /x(AC) = fi{Tnu), 
li{CTnu) = fJ-iC), n{TnuC) = At(T„„), which allow the 
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(16) 



with {^j)mj=i = 0. Figure 4 shows tm as a function of 
equally spaced substitutions C —>■ T„„ and C A. The 
linear scaling of tm{C Tnu) is due to the linear scaling 
of /Uj(C — > Tnu)- The saturating behavior of tm{C — > 
A) is a combination of two opposing factors: the linear 
scaling of ^j, j- (C — > A) and the nonmonotonic behavior of 
/ij (C — > A). Our model also predicts that for sufficiently 
large I 



tm[{CA)i] > tmiQAi). 



(17) 
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This feature is explained by noticing that 
2tm[{CA)i] > tm{C2i) +tm{A2i), See figure 4, which fol- 
lows from the expression for Hj, Eq.(16), and then using 
2t„i{CnAn) « tra{C2i) + <m(^2/), which foUows from 
the linear scaling of tm with iV, for > d in addition 
to Eq.(15). The above findings regarding the behav- 
ior of tm for hetro-ssDNA, fit the experimental results [7] . 

4. CONCLUSIONS 

In the presented model, the translocation of ssDNA 
through a-hemolysin channel is controlled, in addition to 
the voltage, by the interaction between the polymer and 
the pore (A£^j and ^p), and between nearest monomers 
inside the TPP (fj,), as well as by an entropic factors 
originate from polymer segments outside and within the 
TPP. Based on the model, we showed that F{t) can be 
mono or double peaked depending on x and V/Vc- We 
calculated the MFPT to be: r ~ N/{V - Vc), for N > d 
and V/Vc > 1, and tm ~ N/v^iy) for N > d, where 
VxiV) changes from a linear to a quadratic function of 
V with X. In addition we estimated that » lO'^^b, 
and by constructing ^p and for hetro-ssDNA explained 
experimental results regarding the various behaviors of 
tm for hetro-ssDNA. 

An extended version of this model that describes 
translocation through a fluctuating channel structure, 
can be used to describe unbiased translocation, which 
displays long escape times [16]. Translocation of other 
polymers through proteins channels can be described us- 
ing the same framework by changing /i, ^p and AEj. 

Acknowledgements: we acknowledge fruitful discus- 
sions with Amit Meller and with Ralf Metzler and the 
support of the US-Israel Binational Science Founda- 
tion and the Tel Aviv University Nanotechnology Center. 



FIG. 4: tm as a function of the percentage of equally spaced 
monomers substitutions C A and C Tnu, for A'' = 60, 
X = N, z X 1/2, V/Vc = 1.63 and T = 292°K. The curves 
emphasize the effects of the rigidity and the friction on the 
translocation dynamics. See text for discussion. 



APPENDIX A 

We wish to calculate the electrostatic energy differ- 
ence between sates, Ai?J. There are three cases dur- 
ing translocation, which are described by Eqs.(l.l)-(1.3). 
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For any polymer length, Eq.( 1.1) (exit states) describes a 
case for which Nj^cis=0 and Eq.(1.3) (entrance states) 
describes a case for which Nj^trans=^- For the case 
described by Eq.(1.2) (intermediate states), there are 
monomers on both sides of the membrane for a large 
polymer, or no monomers on both sides of the membrane 
for a short polymer. Starting from Eq.(5) we have for the 
entrance states: 



zqV{mj + 1) 



For the exit states AE'J is composed of two contributions. 
One contribution stems from the passage of a monomer 
with an effective charge of zq through the potential V: 

AEP^ = zqV. (A2) 

The second contribution calculated from Eq.(5) is: 



-ijZqV 



Combining the two contributions we find 



AE^ = 



{d — ij)zqV {nij — l)zqV 



d 



d 



(A3) 



(A4) 



when using rrij = d+1 

For the intermediate states and a large polymer, we 
have only the contribution given by Eq.(A2) (the number 
of monomers within the TPP is constant), which can be 
written as: 



rrijZqV 



(A5) 



when using rrij = d which holds for the intermediate 
states. For a short polymer, we have to consider only 
the contribution given by Eq.(5), which leads again 
to Eq.(A5). Eq.(6) is obtained from adding the above 
contributions. 



APPENDIX B 

For calculating A£^|, we start by writing an expression 
for the entropic energy that consists of two terms: 



(Bl) 



Ej -^ represents the entropy cost of storing rrij monomers 
within the TPP and is a linear function of rrij [10]. Ej 2 



originates from the reduced number of configurations of 
a Gaussian polymer near an impermeable wall, and can 
be approximated by [8]: 



\kBTln{Nj^trans), 

^kBTln{Nj^cis), 



-^1,2 = { \kBTln{Nj^transNj^cis 



j<d; 
), N>j>d; (B2) 
j>N. 



Separating the translocation into three regimes, de- 
scribed by Eqs.(l.l)-(1.3), we find that for the entrance 
states AE'I^ is given by: 



while for the exit states 



(B3) 



(B4) 



(Al) with a proportional constant of o(l) 



For the intermediate states AiJ| ^ = because the 
same number of monomers occupy the TPP between ad- 
jacent states. AEj^ for a short polymer has the same 
values as for a large polymer, when adjusting the condi- 
tions for the three cases (exchanging N and d in Eqs.(l)). 

Computing AE^2 from Eq.(B2) results in ^ksTgj, 
where 

ln{l + l/{N -j)), j<d; 
gj = { ln{l + 1/{N - j))(l - - d)), N>j>d; 
Zn(l - l/(j - d)), j > N. 

(B5) 

For a short polymer, gj is similar to Eq.(B5) for the 
entrance and exit states, but gj — for the intermediate 
states because Nj^cis—Nj,trans~0 for these states. Note 
that \gj\ < 1 for all j. Special care is needed when com- 
puting gj for states that belong to the crossover between 
the three situations. For these states a combination of 
Eqs.(B2) was used. From the above contributions we 
obtain Eq.(7). 



APPENDIX C 

To compute t which is given by 

x — l n 

r = ^(-A-i),,, + ^(-A-i),,., (CI) 

5=1 S=X 

we have to calculate the elements of the general inverse 
Toeplitz matrix {A^^)s^x [17]: 

(-A-)...= ^y'^'rr' g:; .<x, (C2) 



A(A)A(A"+i) k 



where A+/_ = [l-|-/-v'(l - 4p+p_)]/2. Substituting the 
expression for p+ and p- into the expressions for A_|_/_, 
we obtain A_|_/_ = which when used in Eq.(C2) 

results in Eq.(12). The summation of each of the terms 
in Eq.(Cl) is straight forward. The first term yields 



/ = a5^(l-y«) = a(x-i-^) (C3) 



where y =^-1^^ and 



1^ A('r)"+l-^) 



a = 



p%A{p 



A(p)A(p"+i)A;' 



(C4) 
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The second term in Eq.(Cl) is: 



r-"-^ ^ f''^ ^ - (n + 1 - x)], (C5) 



1-2/ 



where 



a = 



A(p)A(p"+i)yf 



(C6) 



Combining / and // and rearranging terms results in 
Eq.(13). 

Rewriting Eq.(13) as 



X (n + l)(l-2/-^) 



fcA(p) A;A(p)(l-j/-") ' 



(C7) 



we find that for V/Vc ^ 1, the second term in Eq.(C7) 

vanishes as Keeping the first term in the expression 

for r, Eq.(C7), and expanding to first order in V/Vc the 
explicit form of A(p) 



1 1 + e-^/^-^+i 



A(p) 1 - e-v/^o+-^ V/Vc - 1 
Eq.(14) is obtained. 
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